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1 INTRODUCTION 



Hubble's observations (1929) can be approximated as w oc d , 
where u = cz is the velocity (toward or away from us) of the 
galaxy being observed. Apart from a very few nearby galax- 
ies, redshifts are positive indicating that the galaxies are 
receding away from us. The velocity of recession being pro- 
portional to its distance from us is explained by invoking 
the expansion of the Universe. Hq, the constant of propor- 
tionality, is called the Hubble constant and it measures the 
rate of expansion at the present epoch. The Hubble constant 
enters into various cosmological calculations and its impor- 
tance can never be underestimated. It decides the value of 
critical density pc, the amount of matter and energy required 
to make the geometry of the Universe flat. By comparing pc 
to the observed density one can decide the geometry of the 
Universe. Most importantly it sets the age of the Universe 
(£o) and hence, size of the observable universe {Rob ~ etc). 
Due to its importance determining its accurate value is of 
paramount importance. 

Accurate measurement of the Hubble constant can also 
lead to test the cosmological principal (CP hereafter). Ac- 
cording to CP the universe is homogeneous and isotropic 
at any given cosmic epoch. If the cosmological principal is 
valid than one would expect the average value of Hubble 
constant to be same in different regions and in different di- 
rections. This issue has been addressed by various authors. 
We discuss some of the earlier results below. 
Are we living in a bubble? : If the matter distribution 
is not homogeneous, it causes variation in the value of Hub- 
ble constant. Since gravity pulls we expect that if a region 
of space has higher density then the average density the ex- 
pansion rate will be negatively affected and thus the Hubble 
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constant will be smaller in this region. In contrast to the lo- 
cal mass concentration a regi on with low mass d ensity will 
produce larger value of Hq. (jZehavi et aZ.|[l993 ) provided 
the first evidence for a large local void. They measured Hq 
within and outside 70 /i~^Mpc using SNe la to find that the 
value of Ho was 6.5 % higher than that outside. This indi- 
cates a low density inner region compared to the outside one 
and is known as local bubble or Hubble bubble. The above 
authors had assumed a fiat FRW universe with JIm = 1 in 
their analysis. The variation in the inside and outside values 
of Ho decreases to 4.5% in the ACDM model {Q,m = 0.3, 
S7a = 0.7), however, it does not disappear completely. Re- 
cently (|jha et a/.ll2007l ) revisited the problem using the lat- 
est SNe la data set and detected of the local Hubble bub- 
ble at 3(7 confidenc e level. However in a later publication 
IConlev et alllioOTl ) claimed that it was a misinterpretation 
of color excess of supernovae. At this juncture it is difficult 
to say if the evidence for the local bubble is conclus ive. 
Vari ation in Hq, from HST key data : (McClurc fc Dveij 
l2007h used HST Key Project data (see § [5} to calculate 
the variation in Ho value. The authors find that a statis- 
tically significant variation in Ho of 9kms~^Mpc~^ ex- 
ists in HST Key Project data. The approximate direc- 
tional uncertainty is 10° to 20°. Their results indicate two 
sets of extrema that dominate on different distance scales. 
They find differences as great as ~ 35kms~^ Mpc~^ within 
and ~ 20kms~^Mpc~^ beyond our super-cluster. Within 
70Mpc their results show a statistically significant differ- 
ence of ~ 19kms^^ Mpc~^. This variation does not appear 
to be an artifact of Galactic dust, since there is no consistent 
difference looking in or out of the plane of the Galaxy. In 
fact, the overall structure in the map is inconsistent with the 
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distri bution of dust in the COBE dust maps l|Schlegel et al\ 
Il998l ). 

At this point one pertinent question is to ask "are these 
variations in the measurement of Ho due to a real depar- 
ture from the cosmological principal?" On the contrary it 
is also possible that the data itself has some systematic er- 
rors due to some non corrected physical processes in the 
universe or there could be some real issues with the data re- 
duction/calibration process. In order to comment on "what 
is the real cause of the variations?" , a critical review of the 
measurement methods is required. Measuring accurate value 
of the Ifubble constant is a challenging task. One requires 
accurate measurement of redshift z and the distance d. Al- 
though redshift can be measured with good accuracy from 
the spectrum of the light emitted by the object, the distance 
measurement is difficult. Various methods are employed to 
measure distances, namely TuUy-Fisher relation (TFR), Sur- 
face brightness (SB) fluctuations. Fundamental plane (FP) 
relation, SNe type II, SNe type la, Sunyaev-Zeldovich ef- 
fect (SZE), Gravitational lensing etc. Unfortunately most 
of these methods suffer from the systematic effects arising 
from many different causes. For instance SZE requires the 
3-d distribution/shape of the plasma (hot gas) in the galaxy 
clusters. Radio and x-ray images of the clusters provide only 
the projected x-ray surface brightness and CMB decrement. 
Hence simplified assumptions about the shape of the clus- 
ter are made. Again the assumption abou t the phase an d 
the temperature of the plasma are ad hoc (|Sulkanenll 19991 ). 
In another example the uncertainty in the physical basis of 
TuUy-Fisher relation can cause subtle systematic variances 
in the TF relation with environment. A critical review on 
the subject of Hp measurement methods can be found in 
iJacobv et al] (|l992l '). 

Since there are many sources of systematics, special at- 
tention was needed to measure the accurate value of -ffo- The 
most accurate experiment to achie ve this was the Hubbl e 
Space Telescope (HST) key project l|Freedman et oZ.ll200l[ ). 
We shall discuss the HST key project in the next Section, 
however the issues with this data have been mentioned in 
previous paragraphs. In the present paper we intend to 
put constraints on the CP using the Ho data. This can be 
achieved by looking for direction dependent signatures in 
the data. As has been pointed out earlier detecting the di- 
rection dependent signatures does not guaranty departure 
from isotropy and hence CP. In that case one can constrain 
the r e liabili ty of the data. We use a technique (|Gupta et al\ 
l2008l . I2OIOI ). based on extreme value theory to accomplish 
this. Another important issue with any data set is the pres- 
ence of non-Gaussian errors. Since Central Limit Theorem 
(CLT) predicts that the errors in the data should be Gaus- 
sian, the non-Gaussian errors are undesired and may indi- 
cate some unresolved issues. As a by product of our method 
we are able to detect the presence of non-Gaussian features 
in the errors, which makes our methods useful not only in 
the case of Hq data but for any data set in general. 

Plan of this paper is as follows. We discuss HST key 
project in § [2I which is the main source of our data. We 
discuss our methods in § |4l Since our techniques are based 
on Extreme Value Theory (EVT hereafter) we discuss it 
briefly in § [S] Results and conclusions are presented in § [5] 
and § [6] respectively. 



2 HST KEY PROJECT AND THE DATA SET 

Measuring an accurate value of Ho was one of the moti- 
vating reasons for building the NASA/ESA Hubble Space 
Telescope (HST). Measurement of Ho with the goal of 10% 
accuracy was desi gnated as one of t hree "Key Projects" o f 
the HST ()Aaronson fc Mouldlll98^ : iKennicutt et al.1 [19951 '). 
The overall goal of the HST Key Project was to measure HO 
based on a Cepheid calibration of a number of independent, 
secondary distance determination methods. Many times the 
systematic errors dominate the accuracy of distance mea- 
surement. To overcome this the HST team averaged over the 
systematics and used a number of different methods to mea- 
sure distances instead of relying on a single method alone. 

Determining Ho accurately requires the measurement 
of distances far enough away so that both the small- and 
large-scale motions of galaxies become small compared to 
the overall Hubble expansion. To extend the distance scale 
beyond the range of the Cepheids, a number of methods that 
provide relative distances were chosen. The HST Cepheid 
distances were used to provide an absolute distance scale for 
these otherwise independent methods, including the Type 
la supernovae, the TuUy-Fisher relation, the fundamental 
plane for elliptical galaxies, surface brightness fluctuations, 
and Type II supernova e. The flnal result of HST key project 
jFreedman et al\\200l\ ) was Ho = 72 ± 8km/s/Mpc. 
Data Set : We have c hosen data from the HST key project 
(|Freedman et al\ I2OOII ) as our primary data set. This set 
contains 74 data points and provides a reasonably full sky 
coverage. Different methods used in order to this data set are 
: The TF relation, the FP relation, the SB fluctuations and 
SNe type la, SNe type II. In ad dition we have chosen 2 data 
points from (|Sakai et aZ.ll2000l '). In all the cases recessional 
velocities have been corrected to the CMB frame and thus 
all the Ho values belon g to CMB frame. The fu ll data set is 
published in table 1 of iMcClure fc DveillioOTi ). 



3 EXTREME VALUE THEORY 

In the present paper we investigate the direction dependent 
systematic effects, which exhibit anisotropy in the cosmolog- 
ical data. We identify the direction where the effect of the 
systematics is maximum. To estimate its statistical signifi- 
cance we need to know the distribution of this maximum, 
which can be computed using Extreme Value Theory (EVT). 
Since it is not a common tool in the arsenal of astronomers, 
we begin with a brief introduction of EVT. It was devel- 
oped in parallel to the Central Limit Theory (CLT). While 
CLT describes limiting distribution of partial sums, EVT de- 
scribes how the distribution of extremes (maxima/minima) 
looks like. Below we discuss the theory of maxima, however, 
the results obtained can be easily reformulated to obtain the 
distribution of minima. 

We shall outline the basic ingredients to obtain the the- 
oretical distribution. Let F be a distribution with its right 
end point x* which may be infinite, i.e. 

X* = sup{x : F{x) ^ 1} . 

We randomly choose a sample (Xi, X2, ■ ' ' j^n) of size n 
from this distribution. The maximum of this sample will 
approach x* for large n i.e. max(Xi, ^2, • ■ • ,Xn) — ^ x* as 
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n — > oo. The distribution of maxima has the following prob- 
ability 

P(max(Xi,X2,--- < x) 

= P{Xi X, Xa X • ■ • , X„ < x) 
= P(Xi < x)P{X2 < x) • ■ • P(X„ < x) 
= F"{x) ; 

where Xi, X2, ■ ■ ■ , Xn are assumed to be independent. This 
converges to zero for x < x* and to unity for x ^ x* , which 
means that (x) is a degenerate distribution. A normaliza- 
tion of variable max(Xi, X2, ■ • • , X„) is required in order to 
get a no n-degenerate distr ibution. We choose linear normal- 
ization l|Haan et ffl/.ll2006l ). Let us assume that there exists 
a sequence an > 0, and 6„ real such that 



Gumbel Distribution 



max{Xi,X2 



, Xn ) ^ bn 



has a non-degenerate limiting distribution as n — > oo, i.e 

lim F"{anx + bn) = G{x) , 



(1) 



where ^ G{x) ^ 1, G{x) is the required distribution. It 
is very difhcult to derive the general form of G{x), we men- 
tion some of the results available i n literature. We state a 
theorem due to Fisher and Tippett (|Fisher fc Tippetiri928l ') 
which describes the required distribution. 

Theorem 1. The distribution of maxima G{x : e) has the 
following form : 



G{x\ e) = exp ( — [1 + ex 



-l/s 



(2) 



with 1 -f ex > 0. Where eeR is called the shape parameter 
or the extreme value index. 



Proof of this theorem (| Fisher fc TippetHToii ') is beyond the 
scope of this paper. We only outline a few interesting facts : 



• When e = : 
Tailor expansion of [1 -I- ex] 

G(x) = 



gives e"''. Thus Eq[2]gives 
exp(— e~^) . 



This is known as the Gumbel distribution or extreme value 
distribution of type I. The right end point of this distribution 
is infinity. Also 1 — G(x) ~ as x — >■ oo, indicating that 
the distribution has a thin tail. It is clear from the form of 
Gumbel that it is unbounded on either side that is, there is 
no X for which G(x) = 0. A sketch of Gumbel distribution 
is shown in Fig[l] 

• When e > : 

G(x; e) < 1 for all x, i.e. the right end point of the distri- 
bution is infinity. Also as x — ^ oo , 1 — G(x, e) ~ (ex)"^''^. 
Which indicates that the distribution has a heavy tail. This 
is called Frechet distribution, or type II extreme value dis- 
tribution. One can clearly see that this distribution has a 
lower bound. 

• When e < : 

The right end point of the distribution is — 1/e. This is 
known as WeibuU distribution, or extreme value distribu- 
tion of type III. 

In the above discussion we have not mentioned the loca- 
tion and scale parameters for simplicity. When we introduce 




Figure 1. A sketch of Gumbel distribution. Here we have plotted 
the probability distribution function which is derivative of G{x). 



these parameters, the form of G(x) becomes complicated as 
shown below 



G(x; m, s, e) = exp ( ^ [l + ^{ 



(3) 



where m = 6„ £ i? is the location parameter and s = a„ > 
is the scale parameter. 

For our analysis we shall use the Gumbel distribution 
since the variable of interest there does not have bound on 
either side. Thus our emphasis will be on the Gumbel dis- 
tribution in the rest of this paper, which with scale and 
location parameters takes the form 



G(x; m, s) — exp ^— |^exp — (^^^ — ^ 



(4) 



The probability distribution (pdf) can be derived by differ- 
entiating Eq|4l 



4 METHODOLOGY 



Two new techniques, namely A^2 and Ay statistics based 
on Ex treme value theory, were derived in lGupta et al. 1 (|2008l . 
|2010D . We apply the same method here, however we briefly 
mention these techniques for completeness. 
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4.1 



A^2 statistic 



First we calculate the best fit value of Ho for the complete 

2 



data set by minimizing x' 



defined below 
Hoi — Ho 



(5) 



where Hoi is the i point in the data set and (t; is the 
observed standard error. This minimization gives us the best 
fit value Hot- 

Using this best fit value we now define Xi for the i"' 
data point as 

Hoi — Hob 



Xi = 



(6) 



where Hob is the best fit value of the Hubble constant. We 
shall consider subsets of the data set to construct our statis- 
tic. We define the reduced in terms of Xi as follows 

1 



X 



iVs, 



■E 



X^ 



(7) 



where it should be noted that by 'reduced' we do not mean 
'per degree of freedom', since we do not fit the model sepa- 
rately to the subsets of the data. Here x^ is an indicator of 
the statistical scatter of the subset from the best fit value 

Hob- 

If the Cosmological Principle holds then the value of 
the Hubble constant should not depend upon the direction 
in which it is measured. We use this fact to choose spe- 
cific subsets of data. We divide the complete data into two 
hemispheres, labeled by the direction vector rii, and take 
the difference of the reduced computed for the two hemi- 
spheres separately to obtain Axl, = xLrth - xLuth, where 
we have defined 'north' as that hemisphere towards which 
the direction vector lii points. We are only interested in the 
magnitude of this difference, therefore, we take the absolute 
value of Axrt , and then vary the direction h across the sky 
to obtain the maximum absolute difference 



A^2 = max{lAxrtJ} 



(8) 



We note that since the same data point appears in several 
subsets, the maximization is not done over statistically inde- 
pendent measures of our statistic. Another noteworthy fact 
is that if the direction dependence has a forward-backward 
symmetry then this statistic will not be able to detect it. 
However, due to its ease of construction and use we con- 
sider this simplest of possible statistics. 

To interpret our results we need to know the range 
of A^2 that we can expect if there were no direction de- 
pendence in data, and the noise in the measurements were 
Gaussian. The spatial distribution of measurements is not 
uniform on the sky, therefore, the number of measurements 
in the two hemispheres, for a given direction, varies with the 
direction n in a complicated manner. Therefore one might 
expect the probability distribution function P(A^2) to be 
extremely complicated, however, extreme value theory, § O 
shows that the distribution is, in fact, a simple, two parame- 
ter Gumbel distribution, characteristic of extreme value dis- 
tribution type I: 

A^2 - 1 



■ exp 



exp 



' exp 



A. 





gaussian _ _ 






1 1 





























^2 



(9) 



Figure 2. Here we plot the probability density of A^2 for the 
simulated data. Solid curve represents the bootstrap distribution 
while broken curve represents the theoretical distribution assum- 
ing Gaussian errors. 



where the position parameter m and the scale parameter s 
completely determine the distribution. 

To quantify departures from isotropy we need to know 
the theoretical distribution Pthcory(A^2). Even though we 
know what to expect in a general manner, it is difficult to 
obtain the parameters s and m analytically, therefore, we 
calculate this distribution numerically by simulating several 
sets of Gaussian distributed Xi on the measurement posi- 
tions and obtaining A^2 from each realization. 

If the noise in the data is Gaussian then the above dis- 
tribution adequately quantifies the directional dependence 
in the data. But if the data has non-Gaussian noise then 
the theoretical distribution cannot be used to quantify the 
level of significance of our possible discovery of anisotropy. 
We construct an independent test for directional dependence 
by obtaining the bootstrap distribution Pbs(A^2), which is 
constructed in the following manner. The observed Xi's are 
assumed to be drawn from some unknown, direction depen- 
dent probability distribution. We shuffle the data values Ho^ 
and ai over the measurement positions, thus destroying any 
directional alignment they might have had due to anisotropy. 
Thus we are able to generate several realizations of data and 
estimate the distribution Pbs(A^2). 

In order to know what to expect from this statistic we 
simulate data and calculate the above mentioned Bootstrap 
and numerical distrib utions. These are s hown in Fig. (2] As 
had been discussed in (|Gupta et oZ.ll2010l ) there exists a spe- 
cific bias between the two distributions. Since the numerical 
distribution is obtained by assuming XiS to be Gaussian ran- 
dom variates with a zero mean and unit variance, therefore, 
it does not have any bounds. However the bootstrap distri- 
bution is obtained by shuffling through a specific realization 
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of Xi where the XiS are obviously bounded. It is clear that 
on the average this should produce slightly smaller values 
of A^2 in comparison to what one expects from a Gaussian 
distributed XiS. 

Our results for A^2 statistic in this paper should be 
interpreted with respect to Figure[21 Concerns regarding the 
small number of data points in data set and its effect on the 
efficac y of our method can be addressed (as in I Gupta et al\ 
l|2010ll l by noting that this figure is produced with only 76 
points and the theoretical and the bootstrap distributions 
look similar. 

4.2 statistic 

As mentioned above, Xi does not contain information about 
whether the measurement is above or below the fit i.e. 
greater or smaller than Hot- An obvious generalization that 
does contain this information can be obtained by consider- 
ing a statistic based on XiS. We consider two subsets of data 
defined by two hemispheres labeled by the direction vector 
•h, containing A'north and A'south data points, where the total 
number of data points. A'' = Anorth + A'south, and define the 
quantity 

Ax. = ^ E ^- E ^ • (10) 

Clearly (Axn) = and ((Axn)^) = 1. From the central 
limit theorem (|Kendall fc Stuartlll977l ) it follows that for 
A' ^ 1, the quantity Ax follows a Gaussian distribution 
with a zero mean and unit variance. As in the previous case 
we maximize this quantity by varying the direction h across 
the sky to obtain the maximum absolute difference 

A;, =max{|Axni} . (11) 

This statistic differs from the previous one in that the 
A;^ statistic has a theoretical limit where the position and 
the shape parameters can be determined analytically. Given 
Nd independent directions on the sky we are essentially de- 
termining the maximum of a sample of size A^ where the 
individual numbers are drawn from a Gaussian distribution 
with a zero mean and uni t variance. In the limit Nd ^ 1 the 
parameters are given by (|Haan et fflZ.ll200^ ) 

m = log A^d - log log Nd - log 4tv (12) 

s = 1 (13) 
m 

where we have to additionally assume that the number of 
measurements A' ^ 1, since the distribution for x becomes 
Gaussian only in this limit. This is convenient since at least 
for large data sets, which will be available in the future, 
a comparison with theory becomes simpler. However, for a 
smaller number of measurements (data points) there is a 
possibility that not all directions are independent, in fact, 
it is quite possible that two directions contain exactly same 
subsets in the two hemisphere. In this situation is is clear 
that the total independent directions is a smaller number 
than Nd and thus theoretical distribution would be right- 
ward shifted and also more sharply peaked. For this reason 
we also calculate the bootstrap distribution and the theo- 
retical distribution in the same manner as for the previous 
statistic. 



Table 1. Best fit value for Ho 



Best fit ^ r 

^ -^per doi 

72.0 194.1 2.6 



Theoretical 

Bootstrap 




12 3 4 



a 2 

X 

Figure 3. Here we plot the probability density of A^2. Solid 
curve represents the bootstrap distribution while broken curve 
represents the theoretical distribution assuming Gaussian errors. 
The two curves do not match. for the data lies close to the 

peak of bootstrap distribution. 

5 RESULTS 

First we obtain the best fit value, Hqi,, of the Hubble con- 
stant for the data, which is shown in Table[T] The large value 
of x^ suggests that the error bars may have been underesti- 
mated. 

A^2 statistic : We have applied the A^2 statistic to 
this data and calculated the bootstrap and theoretical distri- 
butions for the HST Key Project data, as discussed earlier. 
The theoretical and bootstrap distributions for this data are 
shown in Fig [3] To interpret our results we compare with 
Fig (2] We see that the Theoretical distribution has been 
shifted far away on the left side of the bootstrap distribu- 
tion. This contradicts to the fact that bootstrap distribution 
should lie slightly to the left of the theoretical distribution, 
as explained above, indicating the non-Gaussian nature of 
errors in the data. HST Key data lies outside the theoret- 
ical distribution, but is close to the peak of the bootstrap 
distribution. This suggests that either the data is free from 
direction dependent systematics or a more sensitive tech- 
nique is required to put constraints on these systematics. 

A;,^ statistic : 

We apply A^,. statistic to the data and calculate the 
bootstrap and numerical distributions. These are shown in 
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Figure 4. Here we plot the probability density of A^. Solid curve 
represents the bootstrap distribution while broken curve repre- 
sents the numerically calculated distribution assuming Gaussian 
errors. The two curves do not match very well. The analytically 
calculated distribution is shown by dotted lines. A;^ for the HST 
key data lies outside one a region around the peak of bootstrap 
distribution. 

Fig [4] Here also the numerical distribution lies to the left 
of the bootstrap distribution, which is in contradiction to 
the expected relative positions of the two. Interestingly, in 
this case we have an advantage, since we can calculate the 
theoretical distribution analytically. This analytic distribu- 
tion is shown by dotted line in Fig [H It was discussed in 
I4.2l that if all the directions are not independent then the 
analytic distribution should lie on the right side of the nu- 
merical distribution and should be peaked sharply. This is 
what we observe in Fig|4l Bootstrap distribution lies to the 
right of the analytic distribution, however, we find that it 
is wider than what we would expect from a true Gumbel 
distribution. This indicates that the bootstrap distribution 
is not truly Gumbel, therefore, it should be compared with 
the numerical distribution, which is calculated in a manner 
that is identical to the bootstrap one. Comparison shows 
that the two distributions differ in a manner identical to the 
difference seen in FigO based on the A^2 statistic; indicat- 
ing a similar level of non-Gaussianity. We also find that the 
position of the data lies outside the 1 a region from the peak 
of the bootstrap distribution. 



dependence in the data, however, statistic which is more 
sensitive to the direction dependence, shows the presence of 
direction dependent systematics at around one a level. 

Acknowledgements Shashikant thanks Arnab Rai 
Chodhuri for supporting part of this work by his DST grant 
(DST0815). 



REFERENCES 

Aaronson, M., & Mould, J. 1986, Astroph. J. , 303, 1. 
Conley, A., Carlberg, R. G., Guy, J., Howell, D. A., Jha, 

S., Riess, A. G., & Sullivan, M. 2007, Astroph. J. Lett. , 

664, LIS. 

Fisher, R. A. and Tippett, L.H.C. , Proc. Cambridge Phi- 
los. Soc. 1928, 24, 180. 

Freedman, W. L., Madore, B. F., Gibson, B. K., et al, 
2001. ApJ, 553, 47. 

Gupta, Shashikant, Saini, Tarun, Deep, Laskar, Tanmoy, 
2008 Mon. Not. Roy. Ast. Soc. 388, 242. 

Gupta, Shashikant, Saini, Tarun, Deep, 2010 Mon. Not. 
Roy. Ast. Soc. 407, 651. 

Haan L. and Ferreira A., "Extreme Value Theory: An In- 
troduction" Springer 2006. 

Jacoby, G. H., et al. 1992, Pub. Astro. Soc. Pacific , 104, 
599 

Kendall, M., & Stuart, A. 1977, London: Griffin, 1977, 4th 
ed., 

Jha, S., Riess, A., G. and Kirshner, R., P., 2007, Astroph. 
J. , 659, 122. 

Kennicutt, R. C, Jr., Freedman, W. L., & Mould, J. R. 

1995, Astron. J. , 110, 1476. 
McClure, M. L., & Dyer, C. C. 2007, New Astronomy , 

12, 533. 

Sakai, S., Mould, J. R., Hughes, S. M. G., et al., 2000. ApJ, 
529, 698. 

Schlegel, D. J., Finkbeiner, D. P., Davis, M., 1998. ApJ, 
500, 525. 

Sulkanen, M. E. 1999, Astroph. J. , 522, 59. 
Zehavi I., Riess A. G., Kirshner R. P., Dekel, A., 1998, 
Astroph. J. 503, 483. 



6 CONCLUSIONS 

We have applied the A^2 and A^ statistics to the HST Key 
Project data. We find that in both the cases the bootstrap 
and the theoretical distributions are very different from each 
other. Thus we conclude that the nature of the errors in the 
data is non-Gaussian. A^2 statistic does not show direction 



